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Abstract. Quarkonia were predicted to be suppressed in the "hot" deconfined matter 
known as the Quark Gluon Plasma (QGP), but they were also seen to suffer from "cold" 
nuclear matter effects (parton shadowing, nuclear absorption...). Both at SPS and 
RHIC, suppression beyond nuclear effects was observed, but the rapidity dependence 
of the RHIC result is not easy to interpret. I review here the current status of these 
i <~"| results, their possible interpretations and the new measurements that could provide 

insights on quarkonia suppression. Some of them were presented at this conference. 

1. A normal introduction 

m 
> 

' Quarkonia suppression, as it was introduced by Matsui and Satz in 1986 [1], should be 

O 

the lighter of which being the J/ip particle, should melt in the quark gluon plasma 



thanks to the screening of the colour charge. Moreover, various bound states having 
different binding energies, they should melt at various temperatures. Measuring different 
quarkonia (J /if), if)' , Xc, the T and Xb families) could thus serve as a thermometer of the 
produced medium, each quarkonia suddenly disappearing above its proper dissociation 
temperature. However, it was soon realized that if the hot quark gluon plasma could 
indeed melt quarkonia, cold normal nuclear matter can also absorbe them. 

The left part of figure [1] summarizes the J/ if) yields (normalized by the Drell-Yan 
process which is known and observed to scale with binary collisions) measured at the 
CERN SPS. From p-p, to various p-A, to S-U, and up to peripheral In-In or Pb-Pb 
collisions, one can describe the data by a simple model in which the J/if> is absorbed 
by the forthcoming nucleons (their number being characterized through the nuclear 
thickness parameter L), with a corresponding cross-section of 4.18 ± 0.35 mb [2]. Both 
the In-In and Pb-Pb more central collisions exhibit further anomalous suppression. 
Moreover, the precise measurements performed by the NA60 collaboration on In-In 
data show that the anomalous suppression saturates in the most central collisions [3]. 
Furthermore, the amount of anomalous suppression is compatible with the fraction of 
J /if) coming from excited states (30-40 %, from if)' and Xc)- One can thus think that, 
at SPS energies, namely around v / s/v r /v — 20 GeV, charmonia behave exactly like the 
predicted quark gluon plasma signature: a fast meltdown of some quarkonia, after a 
simple absorption cross section has been taken into account. 
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Figure 1. Left: J/tp yields normalized by Drell-Yan, as a function of the nuclear 
thickness L, as measured at the SPS [4]. Right: J/ip nuclear modification factors of 
the hottest SPS (Pb-Pb) and RHIC (Au-Au) collisions, as a function of the number of 
participants N part . 



2. The RHIC anomalies 



At RHIC energies, namely a/J/vat — 200 GeV, measurements of J / ip suppression in Au- 
Au collisions by the PHENIX experiment [5] brought up two surprises, as it is shown on 
the right part of figure [TJ First, at midrapidity (red circles), the amount of suppression 
is surprisingly similar to the one observed at SPS (black crosses) if plotted as a function 
of the number of participants N part . There is no fundamental reason for this to happen, 
since the energy density should be higher at RHIC, and the cold nuclear effects could be 
drastically different (in particular, the initial gluons should be further shadowed, having 
smaller Bjorken x). Even more surprising is the fact that, at forward rapidity (blue 
squares), J/ip are further suppressed (by approximately 40%). This is also confirmed 
by a preliminary analysis of the run 7 data shown at this conference. 

So far, we think of two possible explanations of the RHIC J/tp data. 

• First, J/tp could be indeed more suppressed than at SPS, but then recreated during 
(or at the term of) the hot partonic phase from initially uncorrelated c and c 
quarks, the total number of initial cc pairs being larger than 10 in the most central 
collisions [SI El- A large variety of such "coalescence" or "recombination" models 
exists, the latest flavours of which can be found in [El El [El CEH 121 EH Q30 
Qualitatively, they explain why we see less suppression at midrapidity by the simple 
fact that there is more c and c quarks to recombine there. 

• Second, J/tp could be more suppressed at forward rapidity because of nuclear 
effects. Standard gluon shadowing parametrizations do not tend to produce such 
a behaviour [15J but they are poorly constrained by data and further saturation 
effects are not excluded. Interesting exploratory work to derive J/tp production in 

| For seminal publications see references therein. References [5] and [H] were shown at this conference. 
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the framework of the Color Glass Condensate was shown at this conference |lbj . 
It is noticeable that, if such a cold matter effect is responsible for the different 
suppressions observed at mid and forward rapidity, it is possible that even at RHIC 
energy, only the excited charmonia states (if)' and Xc) are melting. 

At this moment, one cannot exclude one or the other of the above scenarios. To draw 
a conclusion, more information is needed and the following section is a comprehensive 
list of observables, some of them having been shown at this conference, that will help 
moving forward in understanding J / if) anomalies seen at RHIC 

3. Moving forward 

3.1. Open charm 

Currently, the open charm absolute cross-section is very poorly known at RHIC. This 
is unfortunate, since it is a basic input to recombination models, the probability to 
reproduce a quarkonia in the QGP depending quadratically on the number of initially 
produced cc pairs. The topic of open charm was covered in another contribution [T7] 
and will not be discussed further here. It is to be noted that the adjunction of new 
vertex detection capabilities in PHENIX and STAR will probably permit much better 
open charm and beauty measurements at RHIC in the future. 

3.2. Cold nuclear matter effects from d-Au 

The effects of cold nuclear matter on J /if) production are also insufficiently known at 
RHIC. Extrapolations from the run 3 d-Au data set were discussed by the PHENIX 
collaboration in [15] and at this conference [18J . Assuming given shadowing schemes, 
namely EKS98 [19] and NDSG [20], and propagating all the uncertainties from the 
d-Au measurements, one derives absorption cross sections of oeks = mb or 

o~ndsg — 2. 21^5 mb. This illustrates that, at present, the uncertainty on such a 
break-up cross section is by no mean better than a couple of millibarns. Furthermore, 
in a data driven approach inspired by [21], PHENIX shows that, if one prefers not 
to rely on a shadowing scheme and to instead propagate the centrality dependence 
measured in d-Au collisions to Au-Au (through a simple Glauber model), the amount of 
suppression due to cold nuclear matter is not know better than shown on figure [2] by the 
shaded yellow bands. Such large uncertainties allow the anomalous suppression to be 
identical at mid and forward rapidity. For instance, the survival probabilities (beyond 
normal suppression) I derive from figure [2] are 55i 3 g % and 38^2 % f° r the most central 
collisions, at mid and forward rapidity respectively. These numbers illustrate that there 
is an anomalous suppression, at least at forward rapidity. They also underline the crucial 
need for a more precise estimation of the cold nuclear effects. Hopefully, it will come 
from the large d-Au data sample recently produced by RHIC, the recorded luminosity 
being thirty times larger than the one the extrapolations of figure [2] are based upon. 
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Figure 2. J/tp nuclear modification factors measured by PHENIX at mid (left) and 
forward (right) rapidity, compared to cold nuclear matter extrapolations from [15] . 



3.3. Transverse momentum 

Since the yield of recombined quarkonia crucially depends on the (unknown) number of 
initially produced cc pairs, the magnitude of the nuclear modification factor itself is not 
an ideal probe of recombination. Instead, the evolution of other observables with respect 
to centrality could carry the signatures of recombined and/or initially produced J/if>. 
We already saw that narrowing the rapidity profile could be attributed to recombination. 
Flattening the transverse momentum would be another of its manifestation. Indeed, the 
spectrum of J / ip formed from uncorrelated, randomly picked, c and c quarks should be 
softer than the one of initially produced J/ if>. Here again, cold nuclear effects should play 
a significant role, by raising the initial parton px through multiple scattering (Cronin 
effect), or even through shadowing [22]. In case of pure Cronin effect, the average pt 
squared should vary proportionally to the amount of centers the partons can scatter 
upon, that can be characterized by L, the thickness parametei{§]. This behavior was 
seen in various p-A collisions at various energies at SPS. During this meeting, the 
NA60 experiment [23] reported a surprising preliminary result exhibiting a slope change 
between p-A at 158 AGeV and previous results, including A-A collisions at the same 
energy. This needs to be confirmed. The left part of figure [3j shows the current status 
of similar measurements at RHIC. Simple linear fits to all the points, including p-p, d- 
Au, Cu-Cu and Au-Au data, exhibit excellent x 2 probabilities and, within these limited 
statistics, there is no evidence for a deviation from the Cronin effect, if any. The mid 
rapidity slope is even compatible with zero while the forward slope has a significance of 
2.7 standard deviations. Their difference is not significant (1.4 standard deviations) and 
deserves more precise data, which will come soon for d-Au (run 8) and Au-Au (run 7). 

In the same spirit, looking at higher pt may provide hints of the processes at 
play in J /if) production. On the left part of figure [3j PHENIX average p\ were 



§ Indeed, in the limit of instantaneous processes, the average number of subsequent or preceding 
interactions are proportional, and the same L thickness parameter can be considered for initial (Cronin) 
of final (nuclear absorption) effects. 
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Figure 3. Left: J/tp average < p\ > versus L from PHENIX. Right: J/tp nuclear 
modification factor in Cu-Cu collisions from PHENIX (low px) and STAR (high px)- 

quoted for 0-5 GeV/c intervals for consistency, but higher p? was reached in the most 
central collisions. At this conference, the STAR collaboration showed a J/ip signal 
(two standard deviations) in Cu-Cu collisions for pt above 5 GeV/c, as well as a more 
significant peak from p-p collisions |24j . While suffering from large uncertainties, the 
nuclear modification factor derived from these measurements is compatible with one, as 
shown by the red stars on the right part of figure El together with low pr corresponding 
PHENIX measurement (empty squares and circles). Statistically correlated high p? 
measurements derived from combining PHENIX and STAR are also displayed (full 
circles and squares). Before to draw any conclusions from this, one should remember 
that this behaviour was observed with high statistics at SPS^JJ] for Pb-Pb [26] and now 
for In-In collisions [23J. In particular, one should keep in mind that the Cronin effect 
may justify this behaviour and will be soon investigated by both STAR and PHENIX 
with the high statistics d-Au data set now available at RHIC. 

The most interesting result of these high px J / if> from STAR, is a clear away side 
correlation with moderate px hadrons (> 500 MeV/c) seen in p-p collisions [23] that 
should start shading some light on the J / ip production mechanisms. 



3.4- Elliptic flow 

Another way of probing the coalescence hypothesis is to measure J/ip elliptic flow. 
Indeed, if charmonia are produced by coalescence of charm quarks, they will somewhat 
inherit their flow, which we know to be quite large from open heavy flavour 
measurements [27J. At this conference, the PHENIX experiment reported a first 
tentative measurement of J/ip's t>2, for Au-Au centrality of 20 to 60 % [28J. As shown 
on the left part of figure HJ these proof-of-principle measurements do not allow one 
to distinguish between models assuming various regeneration (and thus elliptic flow) 
amount. 

|| Hints can also be guessed in the PHENIX Au-Au [5], Cu-Cu [25] and d-Au [15] data, below 5 GeV/c. 
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Figure 4. J/ip elliptic flow measurements. Left: V2 versus px in the PHENIX 
experiment, compared to various predictions (from [28)). Right: Azimuthal 
distribution of J/ip versus the reaction plane, exhibiting clear anisotropy in the NA60 
experiment (from 23J). 

Statistically more significant is the result that NA60 reported, as seen on the right 
part of figure HI It exhibits v 2 = 11 ± 5 % for In-In collisions of centrality 28 to 83 % 
and J/ip of transverse momentum larger than 1 GeV/c. At this energy and number of 
collisions regime, it is very unlikely to be due to elliptic flow. The number of charm 
quarks is low (of the order of 0.05 per average In-In collisions) and they are unlikely 
to thermalize. If confirmed, this anisotropy could arise from an absorption by the 
surrounding anisotropic nuclear matter. In any case, it makes the interpretation of a 
possible J/ijj azimuthal anisotropy at RHIC as a sign of coalescence less obvious. 



3.5. Charmonia excited states 

The normal and anomalous suppression patterns of ip' and Xc higher quarkonia states 
could help solving the J / ip case. In case nuclear effects explain the rapidity dependence 
of J/ip suppression, the anomalous suppression is moderate and compatible with the 
amount of J/ip coming from excited states decays. If so, the excited states should 
mostly vanish in the QGP. On the opposite, if regeneration is at play, excited states 
should also be regenerated and present after the QGP phase. 

Preliminary ratio of J ftp coming from excited states at RHIC energy were shown 
for the first time at this conference [29]. From p-p collisions, the ip' is measured, while 
only an upper limit at 90 % confidence level is obtained for the \ c : 

if) from tp' = 8.6 ± 2.5 % and ip from Xc < 42 % (90 % ULCL) (1) 

Furthermore, one starts to estimate the total beauty cross-section, which translates 
into a fraction of J ftp coming from B at the level of 4^ % [29]. Similar measurements 
in d-Au and Au-Au collisions should bring insights on the J/ip anomalies. 



3.6. Bottomonia 

Beyond charmonia, small samples of bottomonia are now observed at RHIC. This is 
of particular interest, since bottomonia might be easier to understand. First, with 
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much less than one bb pair per central Au-Au collisions, it would be very surprising 
if regeneration plays any role in the beauty sector at RHIC. Second, having higher 
masses, bottomonia are made up of higher momentum partons and will less suffer from 
shadowing (and may even lie in an antishadowing region). 

From the first T signal shown by the PHENIX experiment three years ago |30j . 
a measurement in A-A collisions was very awaited. At this conference, STAR has 
presented a T signal extracted from Au-Au collisions [31]. Together with a companion 
p-p measurement, it should soon provide a nuclear modification factor. For now, the 
collaboration safely quotes an upper limit of RauAu < 1-3 at a confidence level of 90 %. 

3. 7. Onset in the suppression pattern 

From the seminal paper of Matsui and Satz pQ, the ideal signature of deconfinement 
would be a sudden drop of the J /if) yield, when one reaches the dissociation temperature. 
If such a threshold is probably observed in the NA60 data [3] , it is very tempting to see it 
on the midrapidity PHENIX data (remember the left part of figure[2]). This is premature. 
First, such onsets need to be looked at after cold nuclear matter effects are subtracted, 
and we have seen in section [3721 that they are poorly known at present. Second, such 
a discontinuous behaviour does not appear in the forward rapidity data. Last but not 
least, if onset models do fit the midrapidity data, so does a smooth behaviour. However, 
should an onset finally be seen at RHIC, it will immediately be interpreted in terms of 
physical quantities such as a dissociation temperature (see e.g. [32j at this conference). 

3.8. Quarkonia at the LHC 

Finally, even with more information on the above mentioned observables, the precise 
understanding of J/ if) suppression at RHIC may remain delicate, the interplay between 
cold nuclear matter effects, anomalous suppression and regeneration being too intricate. 
If so, one can wonder: what will be the fate of J/if> at the LHC? The number of cc 
pairs will be much higher: likely larger than 100 in a single most central collision. If 
regeneration is at play, then the J /if} nuclear modification factor could even raise with 
centrality, overcoming the normal suppression and the further shadowed initial partons 
(see [8] at this conference for an example of calculation). If so, J /if) enhancement will be 
a spectacular manifestation of reconfinement (and thus of deconfinement). If not, then 
the J /if) yields at LHC may be as hard to interpret as the ones measured at RHIC. But 
LHC will also broadly open the era of T measurements. The very tightly bound ground 
state T could be unsuppressed and thus serve as a reference^! for the excited states (T' 
and T"). This feature should ease the disentanglement of saturation and anomalous 
suppression in the T's family. 

% This will only be feasible once the fraction of T coming from the higher states, including \b, is 
measured or at least taken into account: it is known to be as large as 50 % at the Tevatron |33j . 
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4. An anomalous conclusion 

We saw that J/ip suppression is not yet understood at RHIC, in particular its rapidity 
dependence. To move forward, a lot of tracks can be followed, the most promising of 
which certainly being, on a short time scale, the analysis of the large d-Au sample 
to refine our knowledge of cold nuclear effects. This being said, one should keep 
in view that, assuming conservative cold nuclear matter approaches, some level of 
anomalous suppression is indeed observed at RHIC, characterizing the matter as hot 
and deconfining. 
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